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SUMMARY 


The 2-1/2 kW Advanced Technology Ion Thruster Program 
documented in this report was performed with the goal of establishing 
the technology capabilities and limitations of the 30 cm ion thruster 
design that was representative of the best available technology at the 
time the work was performed (the 700/800 series Engineering Model 
Thruster Design). The information obtained here has been and should 
continue to be useful in determining future technology requirements 
and for assessing 30 cm thruster capabilities in planning the applica- 
tion of ion thrusters for prime propulsion. Data was obtained to show 
the performance characteristics (efficiencies, thermal data, stability, 
etc. ) for a variety of operational parameters in the 0. 5 A to 2. 75 A 
beam current range. These data show that the maximum thruster 
efficiency that can be obtained at any given power level is not critically 
dependent on parameter selection when the effect of beam divergence 
and doubly-charged ions is accounted for. Although the data obtained 
was for a single thruster, the observed effects and trends in parameter 
variation reported here should have general applicability. Relatively 
stable and repeatable operation was demonstrated for this thruster in 
the 0. 5 A to 3. 0 A range of beam currents without any indication of 
thermal or control problems. 

Double ionization in the discharge chamber was investigated 
primarily as a thrust loss factor but also because the doubly- char ged 
ions are thought to be responsible for discharge chamber erosion. 
Discharge parameters such as discharge voltage, discharge power and 
propellant flow split (magnetic baffle current) were varied with the net 
result being that there appears to be a fundamental relation between 
propellant utilization efficiency and the fraction of doubly- charged ions 
measured in the beam for a given thruster (dimensional specification). 

A relatively substantial body of data were collected in support of this 
supposition. The dimensional parameter that most affects the 
relationship between double ionization and propellant utilization is 
found to be the ion optics accelerator grid aperture size (which governs 
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the loss of neutral atoms from the discharge chamber). Examination 
of this effect under this program has only begun to explore this param- 
eter as a possibility for reducing double ionization. Double ionization 
is considered the largest obstacle to achieving higher thruster effi- 
ciency and/or longer lifetime based on the results of this investigation. 

Investigation of thruster-power conditioning interactions showed 
that the 700/800 series Engineering Model Thruster differs little from 
the earlier designs insofar as loss of stable control power conditioning 
requirements and interactions are concerned. Discharge current or 
voltage oscillations are observed in the 10 to Z0 kHz range during typi- 
cal operations, but these oscillations cannot be related to any detri- 
mental thruster behavior. Options for reducing power supply require- 
ments and complexity were studied with the net result being that the 
beam (screen) and accelerator power supplies could be combined, and 
the discharge power supply could be controlled to satisfy the preheat 
and low voltage cathode keeper functions, thereby eliminating the 
requirements tor these power supplies. Employing the discharge power 
supply in this manner would also make more power available for pre- 
heat and thereby shorten the thruster start-up time. 

A study was conducted to determine whether generally appli- 
cable thruster gimbal actuators could be identified that would meet 
most application requirements and present some advantage over other 
choices. The results of this study showed that the gimbal require- 
ments vary widely with respect to the application being considered, 
and there is no well defined "optimum 11 choice of actuator that could 
satisfy this wide range of requirements. Consequently, in place of 
attempting an M optimum 1 1 design, a computer analysis was developed to 
aid in determining gimbal actuator requirements to meet specified 
thrust vectoring and thruster array parameters. 
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INTRODUCTION 


The work described in this report was performed under the 
2- l/2 kW Advanced Technology Ion Thruster Program (Contract 
NAS 3-17831) and is an extension of investigations begun under a prior 
program (Contract NAS 3-16949). The objective of the program was to 
establish the performance characteristics and capabilities of the 700/ 
800 series Engineering Model Thruster (EMT) design that was being 
structurally qualified and endurance tested under other concurrent pro- 
grams (Contract NAS 3-15523 and NAS 3-17803) and thereby to estab- 
lish a baseline for future technology. Much of the data generated under 
this program may be obsolete because of design changes recently 
incorporated into the 900 series Engineering Model Thruster, however 
the trends observed are still valid* 

The work was performed under four distinct tasks. The first 
task was that of documenting thruster performance characteristics 
over a range of operating parameters and power Levels. One objective 
was to find a method for selecting operating parameters so that the 
thruster efficiency is maximum for each power level within the operat- 
ing range* A second objective was simply to document the Limits of 
operation and the efficiency parameters necessary for mission analysis. 
The second task was concerned primarily with identifying parameters 
which control double ionization processes with the goal of minimizing 
the production of doubly-charged ions in the discharge. The third task 
was devoted to investigating power processor-thruster interactions 
with the goal of identifying modifications to thruster or power processor 
design that could simplify or improve the interface. The final task was 
to study gimbal requirements for typical thrust vectoring and to design 
a prototype gimbal actuator on which development for future needs 
might be based. Although only this last task is truly independent of 
the work conducted under the other tasks, this report has been organi- 
zed in sections that describe the work pertinent to each of these four 
tasks. 


II. 


/ 


\ 


DOCUMENTATION OF 30 CM THRUSTER PERFORMANCE 

The performance of an ion thruster is characterized by a 
number of operational properties. Total overall thruster efficiency 
has received primary consideration in the development of discharge 
chamber and ion optical system technology and consequently this param- 
eter has been given most attention. Determination of overall efficiency 
requires measurements of all electrical parameters, propellant flow 
rates, and thrust reduction factors necessary to account for contribu- 
tions to the measured beam current of doubly - c barged ions and no n- 
axial ion trajectories which produce less thrust. At the outset, it was 
considered possible that total efficiency for any given power input could 
be maximized by judicious selection of controllable operating param- 
eters such as discharge voltage or power, magnetic baffle current, 
etc* Thus an algorithm for adjusting operating parameters to provide 
maximum thruster efficiency at any desired power level could be gener- 
ated. While obtaining thruster efficiency data, a number of other per- 
formance parameters were monitored. These parameters included tem- 
peratures of cathodes, isolators, and vaporizers, waveforms of the 
time variations in discharge and beam currents and voltages, and 
erosion rates on selected discharge chamber surfaces. The thruster 
was operated over a range of beam currents from 0. 5 A to 2. 75 A, 
representing a variation in input power from 0. 7 to 3. 8 kW. 

A. Experimental Procedures 

1. Thruster Design and Test Equipment 

The evaluation of thruster performance described here 

was carried out with a thruster which is representative of the 700/800 

1 2 

series EMT design * in all important dimensional, electrical, and 
thermal characteristics. This thruster configuration is shown sche- 
matically in Fig. 1, and is identified as Thruster S/N301-B. The 
thruster tested differs from the EM design only in the materials 
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Figure 1. Discharge chamber schematic. 



(stainless steel instead of titanium), fabrication techniques (fasteners 

instead of welds) and structural elements. Testing was performed in a 

cryopumped vacuum facility, which is 9 ft (2*7 m) in diameter and 18 ft 

(5.4 m) long and maintains an ambient pressure in the range from 1 to 

5 x 10"^ Torr in the vicinity of the thruster during testing* The thruster 

3 4 

is operated and controlled using a transistorized power processor ' 
with principal power supplies and metering connected as shown sche- 
matically in Fig* 2* All evaluation was performed with power proces- 
sor control loops closed, as would be normal in an actual flight appli- 
cation* Each set of operating conditions was maintained for a sufficiently 
long time period to obtain propellant flow rate data. 

2, Electrical Measurements 

Electrical measurements were of two types: steady 

state meter readings for power and thrust determination, and oscillo- 
scope recordings of beam and discharge voltage and current waveforms 
for stability evaluations. Time-varying parameters were monitored 
only in the dc to 50 MHz range since higher frequencies are not consid- 
ered important to the operation or control of the thruster system. Alt 
metering is calibrated with respect to primary standards, and metered 
electrical measurements are accurate to within ±1%. 

3. Propellant Flow Measurement 

Propellant flow rate was measured by recording pro- 
pellant reservoir content at regular time intervals and by graphically 
obtaining a "best fit" linear rate to the recorded data. Two methods 
were used to monitor propellant reservoir content: one measures the 

volume of mercury in a calibrated glass pipette while the other weighs 
the mercury stored using load cells calibrated with primary standards. 
These methods were used interchangeably and produced results which 
agreed within 1%. 



Figure 2. Schematic diagram of principal thruster power supplies 
and meter locations for direct measurement of thruster 
performance parameters quoted in Table I. 
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Correction Factors and Beam Probes 


Thrust corrections which account for doubly-charged 
ions and nonaxial beam trajectories were measured with a collimating 
mass spectrometer probe* This apparatus measures the singly- and 
doubly-charged beam current components as a function of angle and 
grid location along a diameter of the ion extraction system. Data are 
recorded automatically and analyzed by computer. Assuming azimuthal 
symmetry, the measured current densities can be integrated and thrust 
factors are computed for use in correcting the measured beam current 
so that the thrust can be calculated from beam current and voltage on 
the basis of b beam of singly-charged paraxially ejected ions. These 
correction factors are defined as follows: 


T = aF,T 
a t c 


(i) 


where T is the thrust calculated from measured beam current, a is the 
c 

factor which accounts for the reduced thrust produced by doubly-charged 

ions, accounts for the thrust reduction on nonaxial or diverging ion 

velocity components, and T is the actual thrust produced. The thrust 

a 

determined by using these measured correction factors with metered 

values of beam current and voltage has been reported to agree to within 

1% of the value determined by direct thrust measurements at NASA 

5 

Lewis Research Center. Another parameter has been defined which 

relates the ionized mass flow rate m, to the neutral mass flow rate m 

l o 

by 


m. I-qM 

• — S • 

m K e m 
o o 


( 2 ) 


where Ig is the beam current, p is the factor which takes into account 
the doubly-charged ions, e is the electronic charge, and M is the mass 
of a mercury ion. In addition to the ExB separator, a collimated 


Faraday cage was also used for the investigation of fringe ions and 
charge-exchange ions. Since this probe could be positioned to much 
larger dispersion angles than the ExB assembly, it allowed a more 
detailed examination of beam divergence characteristics. 

5. Temperature Measurements 

Thermocouples were attached to thruster comp >nents in 
12 locations which have critical temperature requirements. Tempera- 
tures were monitored at regular time intervals, as in the case of the 
propellant flow determination in order to make certain that steady- state 
conditions had been established. It should be emphasized that tempera- 
ture distributions reported here are accurate only for steady state 
thruster operation in a large, cryopumped vacuum facility. 

6. Erosion Measurements 

A numer of techniques were attempted to monitor ero- 
sion of discharge chamber components by ion sputtering. These 
included comparison of component photographs between tests, use of 
small cover plates which had alternating thin layers of different colored 
metals (iron and copper for example), and cover plates which had been 
highly polished and then coated with appropriate materials (e. g. , tanta- 
lum, iron). Only the latter technique using a Dektak profilometer to 
measure the depth of material removed from a coated slice of polished 
silicon wafer produced any quantitative information. 

B. Experimental Results 

Parameter variation followed a schedule whereby baseline per- 
formance conditions were established at each beam current value 
(0* 5 A, 1.0 A, 1. 5 A, and 2, 0 A) for fixed values of the following 
parameters; 


Beam Voltage, V 

1100 

Accel Voltage, V 

500 

Discharge Voltage, V 

37 

Magnetic Baffle Current, A 

3. 5 


22 





Having obtained all forms of data for this set of parameters, a 
systematic variation was pursued to determine whether increases in 
total corrected thruster efficiency could be produced and to document 
which other parameters might be influenced. The more significant 
results of these experiments are described in this section* 

1. Performance Mapping 

Since the foremost objective of this parameter variation 
study was to maximize the total corrected efficiency for thruster opera- 
tion over a beam current range from 0* 5 A to 2* 0 A, this aspect of the 
experimental results will be discussed first. Total corrected thruster 
efficiency, is defined here as 


T ~ 


T^ 

a 2„2 

2m P = n p T lu 0 ' F t 
o r 


(3) 


where 


Vb 


^u m 

o 

and P is the total input power, and Ig are the metered beam voltage 
and current, respectively, and the remaining quantities have been 
defined above* At the outset of this investigation it was considered 
possible that tradeoffs might exist between the quantities in eq* (3) 
which could appreciably affect the value of the total efficiency- The 
performance data which were obtained to establish a basis for selecting 
optimum parameters are presented as Table I. These data were col- 
lected over a number of moderate length testing periods and are 
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representative of operating points which are consistently reproducible* 
Higher and lower efficiencies have been obtained with both the thruster 
used in this study and other thrusters of the EMT design. When all 
sources of error are considered, the accuracy of the overall corrected 
efficiency is limited to ±1%. The data in Table I have been organized 
for comparison of variations at each beam current level. "Normal 11 
neutralizer parameters have been used in the computation of r| , r| , 
and r\rjr because the neutralizer used in these tests was not operating 
with what are considered "representative" characteristics* One of the 
first observations to be made regarding the data in Table I is that for 
any given beam current, there are a number of choices of the controlled 
parameters for which the variation in is not greater than the abso- 
lute accuracy of rj^. Consequently, it appears that rj^ is relatively 
insensitive to the selection of control parameters, and some factor 
other than such as minimizing the ratio of doubly-charged ion cur- 

rent to total current, or simplifying of control algorithms could be used 
to make the parameter selection, 

A number of trends or tendencies can be explored by cross- 
plotting data from Table I, The dependence of on power input, or 
the throttling characteristic is considered first. Figures 3, 4, and 5 
show this "power throttling" characteristic for several combinations of 
parameters. It is apparent, as stated earlier, that the choice of param- 
eters, or an inadvertent variation of control parameters, does not sig- 
nificantly alter this characteristic. It is interesting to compare actual 
thruster data with a theoretical throttling characteristic computed on 
the basis of Kaufman's maximum propellant utilization theory, ^ In 
brief, this theory states that "the loss rate of un-ionized propellant at 
the maximum utilization is nearly a constant over a wide range of total 
propellant flow rate, permitting rapid prediction of total flow rate 
effects from operation at one total flow rate, " Quantitative agreement 

between thruster data and performance predicted by this theory has not 

5 

always been good, possibly because the presence of doubly-charged 
ions has not been adequately taken into account. Similarly, it has also 
been considered that this theory is invalidated by the flexibility in 
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Table I. Engineering Model Thruster Performance Data 
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operation offered by the use of a magnetically controlled baffle which 
provides relatively independent control of discharge voltage, current, 
and cathode propellant flow rate. Since the 3 0 cm thruster was opti- 
mized at the 2 A beam current level, an optimum data point"' has been 
selected from Table I to compute the loss rate of un-ionized propellant, 
by the relation 

= m o - IgP = 0.274 A (equivalent) . (4) 

The theoretical curve is then generated by using the values of q , & , (3, 

P 

F^, and P from Table I for each beam current, 1^, and calculating the 
total corrected efficiency by the equation 

= (pi B + 0.274) \ Q F t * 

Figure 6 shows the remarkably good agreement obtained when the curve 
generated using eq, (5) is compared with measured values for Table I. 
Thus, the phenomenological validity of Kaufman's propellant utilization 
theorem is supported by the results obtained here. If one accepts this 
theorem, then the throttling characteristic generated by eq. (5) repre~ 
sents the maximum achievable characteristic and the variations shown 
in Figs. 3, 4, and 5, are expected to have equal or lower values, at 
best. 

2. Magnetic Baffle Effects 

It can also be seen in Table I, that for any given discharge 

voltage and discharge loss value, the cathode propellant flow rate is 

directly proportional to magnetic baffle current. Choice of magnetic 

baffle current is based on stabilizing discharge operation to prevent 
7 

mode shift (whereby the discharge voltage is no longer controllable 
in the "normal" range by cathode flow rate). Optimum propellant 

4 First entry in Table I at 2 A beam urrent and using 32 mA equivalent 
neutralizer flow (best point from neutralizer mapping). 
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efficiency is consequently obtained at the minimum cathode propellant 
flow rate which is greater than the instability limit and control of the 
magnetic baffle to satisfy this condition ranges from 30 to 40 ampere- 
turns for a 5.4 cm diameter baffle. 

3. High Power Operation 

As noted in Table I there are three data points which 
were obtained at beam current levels above the nominal 2. 0 A. These 
points were considered precursors to the 4. 0 A beam operation, and 
served to identify any obvious thruster problems which may need atten- 
tion. Thus, a little more discussion is in order concerning operation 
at these power levels. 

At both the 2. 5 A and 2. 75 A beam current levels the thruster 
operated satisfactorily with no apparent problems such as excessive 
temperature excursions or shorting of the accelerator grids. Due 
to power conditioning limitations, operation at 3. 0 A beam current 
was maintained for only a short time, but the thruster itself appeared 
to operate quite well* At these higher power levels, cathode tip tem- 
peratures remained in the 900 to 1000°C range and isolator tempera- 
tures were within the 250 to 300°C range. Although the thruster could 
continue to operate at 185 eV/ion, the total extraction voltage had to be 
increased to the 1750 to 19 00 V range in order to avoid operation near 
the perveance limit. Judging from the periods of thruster operation at 
higher power levels there appears to be no significant reason why con- 
tinued operation at these and higher power levels cannot be maintained. 
The amplitude of the discharge oscillations continued to decrease at 
higher beam currents, and as a result both the bea„m and discharge con- 
trol functions were very stable. 

4. Thermal Data 

In the process of collecting the performance data shown 
in Table I, steady- state operating temperatures were recorded at the 
locations shown in Fig. 7. Table II gives typical values for these tem- 
peratures for a number of the operating points given in Table I, 
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Neutralizer thermal data has again been omitted because the neutralizer 
used was not suitably representative of the EMT design. The only 
component temperatures which are not totally satisfactory during typical 
operation are the isolator temperatures at low power operation {<1300 W). 
These components require heater power {normally used only during 
startup) to prevent propellant condensation. Cathode tip and insert tem- 
peratures remain in a satisfactory range even up to the 2. 5 A beam 
current level. It is anticipated that all of the component designs would 
be satisfactory from a thermal standpoint even at higher beam currents* 

5. Oscillations 

One characteristic of thruster performance which cannot 
be adequately explained with regard to either cause or effect is the 
oscillatory behavior observed in the discharge chamber. These oscil- 
lations are manifested primarily as a modulation of the discharge cur- 
* rent with lesser modulation of discharge voltage and beam voltage and 

current. The amplitude and frequency of this modulation depends on 
- several factors such as magnetic baffle control current, cathode pro- 

pellant flow rate, discharge voltage, and configurational parameters 
which affect the magnetic field in the cathode orifice region. Oscilla- 
tion amplitude and frequency can also be influenced by power processor 

g 

output characteristics. Typical waveforms of discharge current for 

3. 5 A magnetic baffle current at four beam current levels are shown 

in Fig. 8. It can be seen that the peak-to-peak modulation amplitude 

remains approximately the same at each beam (discharge) current, 

thereby increasing the percentage of modulation with decreasing beam 

(discharge) current. The fundamental frequency at full beam power 

operation i s approximately 20 kHz and the oscillation is quite periodic 

and stable. Lower frequencies appear and tend to dominate at lower 

> discharge currents. The current is very nearly interrupted at a 2. 5 

kHz rate at the 0. 5 A beam point. Synchornized waveforms of dis- 

y charge current and voltage indicate that these quantities are out of 

phase by approximately 90 to 180 degrees. Very little more can be 

said about these oscillations, although operating conditions have been 

8 

found which essentially eliminate them. No evidence has been 
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established to substantiate any loss of efficiency, or increased 
doubly-charged ion oroduction which can be associated with these 
oscillations. These oscillations are somewhat troublesome to power 
processor and control functions, however, especially at the lower beam 
current operating points. Except at lower beam operation where dis- 
charge current is nearly interrupted, the beam current fluctuations are 
observed to be negligibly small and uncorrelated with discharge cur- 
rent oscillations. This indicates that the observed oscillation is more 
a variation in electron transport process than in plasma density or 
ionization process. The only conclusion that can be reached at this time 
based o'u this observed oscillatory behavior is that extension of thruster 
operation to higher power levels, rather than lower power levels, would 
be preferable. 

6* Erosion Monitors 

Erosion of thruster components was monitored through- 
out the performance mapping tests by several techniques and in a num- 
ber of locations. Since backsputtered deposition from the steel beam 
target is a competing process, only areas where sputter erosion is 
relatively severe were observed to display material loss, most of the 
locations showing either no loss or a coating. Careful analysis of these 
coated surfaces would be required to deduce erosion rates and this was 
not attempted. Consequently, only two thruster components could be 
quantitatively analyzed for erosion rate. The neutralizer cover plate 
was monitored using an n erosion badge" on which eight layers of cop- 
per and chromium were deposited in alternating layers of 1000 A and 
250 K, respectively, (This provides approximately equal sputtering 
times per layer. ) Erosion was observed from the downstream-most 
tip of the neutralizer cover to the point at which the thruster cover 
shields the neutralizer from the ion beam. Erosion rate was equiva- 
lent to 50 A/hour for chromium. The other location which showed 
measurable erosion was the downstream surface of the baffle. Initial 
attempts to monitor erosion of this surface using the alternating layer 
technique were unsuccessful because the layers were completely 
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removed in the test times allotted between inspections. A successful 
semiquantitative measurement was completed by fastening a thin, highly 
polished and coated silicon wafer to the baffle surface. The coating 
on the wafer was a sputter-deposited layer of steel, 6000 A thick. The 
wafer was masked in several areas to provide ,r bench marks’ 1 for obtain- 
ing a depth measurement with a precision profilometer . After 20 hours 
of operation, the step measured from the tench mark region was 
approximately 7000 A, or 350 A eroded per hour. Several repetitions 
of this experiment at different baffle locations yielded erosion rates 
from 200 o 500 A per hour with erosion rates being highest near the 
baffle center and edge. These measurements are considered to be mini- 
mum erosion rates since erosion and deposition are competing processes 
and the silicon wafer is suspected of charging somewhat because of its 
lower conductivity. On the ba^.s of these rates, the baffle should not 
be eroded away until approximately 20,000 hours of operation, which 

is about a factor of three greater than observed in the endurance test 

q 

described by C. Collett. Consequently, it must be concluded that for 
measurement of sputtering erosion using the techniques described 
above, greater care must be taken to prevent charging of the sample 
being studied, and that backsputter ed beam target deposition must 
either be eliminated (through use of a frozen mercury collector) or 
accounted for in some indirect manner. 

C . Additional Results 

1. Neutralizer Studies 

Prompted by abnormalities chat were observed during 
the performance mapping, and in the endurance test of thruster S/N 701, 
further investigation of an 800 series neutralizer was conducted. The 
purpose of this phase of testing was to characterize a new neutralizer 
(S/N 806) and to observe its characteristics with time so that a basis 
for future remedial action could be established. Control curves for 


the new neutralizer are shown in Fig* 9(a) for various beam current 
levels, where keeper voltage is plotted as a function of vaporizer tem- 
perature. Figure 9(b) shows the variation in floating potential for these 
same points* Throughout all tests the neutralizer keeper current was 
maintained at 1,5 A* After approximately 6 hours of operation the 2. 0 A 
beam condition was repeated* After it was found that the neutralizer 
characteristic had shifted to a somewhat higher keeper voltage level, 
operation was then continued at full power to observe any further 
changes in performance* The curves of Fig* 10(a) show the various 
keeper voltage functions after accumulating the time as indicated, while 
Fig. 10{b) again shows the floating potential variations. In general, 
there seems to be a decrease in the slope of the control curve with 
time, and after approximately 32 hours of operation a more dramatic 
shift in voltage was noted. This particular higher voltage character- 
istic was first obtained when high frequency power conditioning was 
substituted for the 60 Hz supplies. However, a return to the low fre- 
quency supplies did not restore the preceding control characteristic 
and the neutralizer continued to operate on this higher voltage control 
characteristic* 

The variation in control characteristics exhibited by this 
neutralizer was very similar to the characteristics observed with 
the life test neutralizer S/N701 where two modes of operation were noted. 
In the life test the shifting appeared to be random and operation in any 
one of the modes would continue for many hours. The less dramatic 
shifting of the keeper voltage versus propellant flow characteristic has 
been noted in all 700 and 800 series neutralizers; in particular, the 
shake test neutralizer (S/N 703) displayed control characteristics simi- 
lar to those of S/N 806 except that no mode requiring a higher keeper 
voltage set point (~16V) was observed. 

Attempting to identify the source of the observed control 
variations one notes that there appears to be a rather steep thermal 
gradient along the cathode tube between cathode tip and insert location. 
Figure 11 shows a schematic drawing of the neutralizer and designates 


















5047-23 


ROLLED TA 



Figure II. Cross section of neutralizer cathode indicating 
location of thermocouples. 
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the thermocouple locations. The temperatures recorded at these 
locations during the performance mapping of the neutralizer are given 
in Table III, It appears likely that insert temperature is too low in the 
desired operating regime to contribute to either emission processes or 
aid in supplying barium to the higher temperature surfaces downstream 
of the insert. These results prompted a more in depth investigation at 
NASA LeRC that resulted in specification of a barium aluminate impreg- 
nated porous tungsten insert to replace the coated foil insert. 

2. Beam Studies 

The evaluation of the effluent ion beam from a 30 cm 
thruster was prompted by a decrease in mass utilization efficiency 
which was noted during the course of the performance mapping. Upon 
returning to 2. 0 A beam current after considerable operation at 1/2 
and 1/4 power levels, it was found that the mass efficiencies had 
decreased 1 to 2%. An examination of the grid set revealed an increase 
in the diameter of the accelerator apertures at increasing distance from 
the thruster centerline. As can be seen from Table IV, the accel aper- 
ture size varies in ar. almost linear fashion with radius. Because of 
the 0, 5% screen grid compensation, the ’Vectoring" at the edge of the 
grid is considerable and results in ion crossover in the aperture region. 
Also, at lower power levels the beam divergence is significantly greater 
than at full pc^'^r. as is shown in Fig. 12 where the normalized 
ion current density is plotted as a function of dispersion angle. Both 
the 0. 5 A beam current and the 2. 0 A beam current curves are repre- 
sentative of a radial location of 10.25 cm. Although the absolute cur- 
rent density is less for the 0. 5 A condition, operation at these low beam 
levels for extended periods could result in increased aperture diameter 
due to direct interception. Based on an average enlarged hole size of 
1. 68 mm, as compared to the 1. 58 mm original size, the equivalent 
neutral loss difference and subsequent drop in mass utilization could be 
calculated. Although the calculation was approximate, the indication 
was that the enlarged accel holes could certainly be responsible for the 
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Figure 12, Single ion beam divergence at r = 10.9 cm for 
beam currents of 2 . 0 and 0.5 A. 






Table IV. Variation of Accel Aperture Diameter 
Across Grid S/N 644 
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Thruster 
Radius* cm 



— 
- 15 

-12.3 

- 10 
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performance drop. Optics S/N 650, which are identical to optics 644 
except that they have uniform aperture size, were substituted and the 
thruster performance was again evaluated. As anticipated, the mass 
utilization as well as the perveance were restored to the original design 
conditions. Table V shows the performance parameters for the two 
optics sets. Based on these results, it appears that the total accelerat- 
ing voltage must be reduced for lower beam currents to preserve the 
focusing of the individual beamlets in order to minimize any erosion of 
the accelerator electrode apertures which result from direct interception. 

In conjunction with the beam profiles done for aperture size 
determination, two more detailed studies of high angle beam trajec- 
tories were also done. These studies included primary as well as 
charge exchange ions. Most of the primary beam ions are well colli- 
mated with more than 99% of the beam current contained in a 20° half 
angle. However, there is a small class of ions, most of which origi- 
nate near the periphery of the grid which are emitted with angles much 
larger than the nominal 15° half angle. To measure these fringe ions a 
narrow angle collimating probe was built to sample the angular current 
density distribution emanating from a small area on the accelerator 
grid. Figure 13 shows the geometric configuration of this collimating 
probe, which allows it to travel in circular arc while viewing the same 
grid area. This probe has a 2° acceptance half angle resulting in a 0,7 
cm diameter viewing area at a probe-accel distance of 1 0 cm. Figures 
14 and 15 show probe characteristics taken at different accel voltages 
at the central and edge regions of the accelerator grid. These curves 
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PROBE TRAVELS IN ARC, 
VIEWS SAME AREA ON GRID 

COLLIMATING PROBE 


Figure 13, Geometric configuration of collimating probe used 
to sample angular current density from accel 
electrode . 







Figure 15. Experimental data showing angular 
current density distribution ema- 
nating from edge holes showing 
fringe beam ions at high angles. 
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Table V. Performance Parameters for Optics S/N 650 and S/N 644 


Optics 

Weighing System 

Pipette 

T) % 
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mam 

rh j i 
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rh 
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S/N 644 

•93.4 

2. 045 

0. 094 

94.2 

2. 03 

0. 092 

S/N 650 

94. 3 

2. 018 

0. 102 

95.2 

2. 00 

0. 100 


T1871 


show how the beamlets from the central region have appreciably less 
angular divergence than the edge beamlets which originate from a more 
concave plasma sheath boundary* The traces also show how the pro- 
files become more peaked as the accel (Y ) voltage is reduced. The 
current trace taken at the edge shows a small group of ions which are 
emitted with angles up to 80°. These ions were not observed until a 
1 cm wide slit was cut in the ground screen electrode. The ground 
screen geometry and placement is at the correct location to intercept 
these ions* Usually after a few hours of operation the edge of the 
ground screen facing the edge row of accelerator holes has the charac- 
teristic clean appearance of a sputtered surface* Retarding field 
measurements on these "fringe” ions showed that they had the full beam 
energy; they were not poorly focused Group 2 charge exchange ions. 

It was also observed that the height of this current peak decreased with 
accel (V ) voltage. First attempts at modeling these fringe ions 
assumed they were crossover trajectories from the corner of the 
plasma sheath. Computer trajectories showed that it was not possible 
to simulate the 80° trajectories even with an ion with a 30 eV initial 
energy directed in the outward radial direction. The fact that the 
fringe ions were not observed when the accel voltage was reduced to 
-400 indicated that these ions might be related to the end effect at the 
beam edge where the virtual ground neutralization plane curves away 
from the accelerator grid. The field lines between the accel and this 
virtual ground boundary are in a direction to bend the peripheral ions 
in an outward direction. As the accel voltage is reduced this force 
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decreases. Figure 16 shows the beam edge geometry used to model 
these trajectories. Trajectories were run from the simulated accel 
electrode making angles of 20° with respect to the edge hole centerline. 
The assumed neutralization plan was placed a distance of 0*05 in. from 
the accel electrode* These calculated trajectories do exhibit the behavior 
observed experimentally. It is seen that the edge trajectories exit with 
angles from 60 to 80° in agreement with the experimental data. 

The other experiment performed with the collimator probe was 
to measure radially directed charge- exchange ions formed in the down- 
stream beam plasma. The probe was mounted perpendicular to the 
beam centerline and could be moved from a position 4 cm upstream 
to a position 15 cm downstream of the thruster ground screen. The 
probe was positioned a distance 6. 5 cm outward from the nominal 15 cm 
beam radius. The experimental measurements in Fig* 17 show that as 
the probe is moved downstream from the ground screen electrode, the 
current rises from zero to a peak value at the axial position of the grid 
centerline. In this region the probe is viewing increasing surface areas 
of the dished grid which extends outward about 1 cm from the ground 
screen. Downstream of this peak the current slowly decreases with 
increasing distance. Comparing the above results with data taken at a 
probe angle of -15° (in order to view ions exiting perpendicular to the 
nominal 15° beam spreading angle) we see that the probe current is 
down more than an order of magnitude* This indicates that the flux of 
the Group 4 charge exchange ion is directed predominantly in the radial 
direction rather than in a direction perpendicular to the beam envelope. 
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Figure 16. Model used to simulate fringe beam 
ions which are deflected outward due 
to the distortion of the virtual ground 
neutralization plane. 
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Figure 17. Measured current of radially directed 
Group 4 charge exchange as a function 
of distance downstream from ground 
screen electrode for 0 to 90° and 
105° 



III. 


DOUBLY- CHARGED ION INVESTIGATION 


This section describes the work which was performed to 
evaluate multiple-ion production and to identify means by which these 
ion concentrations could be reduced. The detrimental effects of the 
doubly-charged ion current are manifested in two ways: 

1* By an erosion of discharge chamber components, 

in particular, the baffle and screen grid 

2. By a reduced thrust per unit current that results 

in lower total efficiency. 

The first effect is now considered of primary importance 
although it was the reduction in efficiency that originally motivated this 
investigation. The intent of this section, therefore, is to describe 
briefly the apparatus used in this phase of testing and then to summarize 
and discuss the experimental results that were obtained* 

A. Beam Probe and Measurement Technique 

The ExB mass separator probe was designed and built under 
previous contracts (NAS 3- 16528 and NAS 3-16949). References 8 and 
9 describe in detail the entire probing system thus only a brief descrip- 
tion of the probing technique and method of data analysis will be given 

here* Because of the 0. 6° probe acceptance angle, the viewing area 

2 

of the probe at the accelerator grid is approximately 0. 132 cm . Pre- 
programmed stepper motors position the probe to view selected radial 
locations in the thruster diameter from various dispersion angles 
ranging from +40° to -40°, Collector currents for both the singly- and 
doubly-charged ions are obtained from seven dispersion angles at each 
of four radial locations with the results stored on paper tape. These 
results provide input data for the computer where the double to single 
ion current is calculated by averaging overall dispersion angles and 
radial locations. In addition, the program provides the following 
information. 
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I. 

Thrust and mass efficiency correction factors 
a, to account for the presence of doubly- 

charged ions* 

h 

2. 

Thrust correction factor to account for beam, 
divergence F t * 


3. 

Plots of the radial profiles of the calculated 
current densities. 



These factors “the data acquisition system, the number of 
data points used, and the associated averaging effects — make this 
technique repeatable to within 2%. 

B. Small Hole Accelerator Grid 

At the outset of the program it was felt that small hole accel- 
erator grid (SHAG) optics would be beneficial in achieving thruster 
operation at reduced double ion concentrations . Because of the long 
lead time necessary for grid manufacture, a set of SHAG specifications 
was a priority item before any testing was to begin, A minimum accel 
aperture diameter was desirable from the performance viewpoint; how- 
ever, the holes should not be so small that considerable time would be 
spent in machining apertures with the added possibility of destroying 
beam divergence properties by incorrect compensation. 

In an effort to establish a minimum accelerator grid open area 
fraction, two techniques were employed* One technique was to cover 
apertures in an existing grid set with thin Ta foil, and permitting the 
beamlet to machine a new aperture* Thruster operation at 2. 0 A beam 
current quickly sputtered the foil and an aperture diameter of 0, 089 cm 
was determined as minimum on the thruster axis. A second technique 
was to compute the ion trajectories for various accel aperture diam- 
eters* Results of the computer analysis are shown in Fig* 18, whereby 
a reduction of diameter from the present 0. 156 cm to 0. 103 cm would 
not result in any direct interception and agrees well with the experi- 
mental data. If such a grid were manufactured the open area fraction 
would be approximately 20%, Table VI summarizes the parameters 
which were chosen for the manufacture of the SHAG optical system. 
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Table VI* Parameters for Low Transmission Accelerator 
Grid Ion Optical System Design 


Screen Grid 

in. 

cm 

Thickness 

0. 015 

0. 038 

Aperture Diameter 

0. 075 

0. 191 

Aperture Spacing 

0. 087 

0.221 

Reduction in Spacing Pattern, % 


0.4 

Open Area, % 


67. 6 

Accelerator Grid 

in. 

cm 

Thickness 

0. 015 

0. 038 

Aperture Diameter 

0. 040 

0. 102 

Aperture Spacing 

0. 037 

0. 221 

Open Area, % 


19.2 


T187Z 


C. Experimental Results 


Two approaches have been empLoyed to reduce doubly-charged 
ion densities while still maintaining high total efficiency; 

1, Parametric variations establish optimum operating 
conditions with minimum double ion content* This 
phase of testing was carried out in conjunction with 
the performance map described earlier. 

2. Thruster geometry changes, in particular the ion 
opticaL system, were evaluated in order to estab- 
lish optimum performance, again with minimum 
double ion content. 

1. Parametric Investigation 

As the performance mapping of Section II has -shown, 
the total thruster efficiency as a function of input power is relatively 
insensitive to the selection of control parameters. Thus, the per- 
centage of doubly-charged ions should be weighed more heavily in the 
selection of these control parameters provided the necessary relation- 
ships are known. The performance mapping provided a range of 
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thruster operating conditions on which double ion studies could be 
based. The range of parameter selections were 


V D IS 35 to 41V 


185 to 250 eV/ion 


I Q 0. 5 to 2. 75 A 

ID 


V NET 1100 to 1300 V 


Selected points from the performance mapping of Task I are 
shown in Fig. 19 where the double ion current density ratio is plotted 
as a function of corrected mass efficiency. The soitd line represents 
the data points which would be obtained by varying the beam current 
at a discharge voltage of 37 V with discharge losses of 185 eV/ion. 

The data appear to have the following trends: 

1. At a given beam current, any attempt to increase 
mass utilization efficiency by increasing the dis- 
charge voltage results in an increase in the 
doubly-charged ion ratio with little or no increase 
in the corrected mass efficiency. The same 
appears to be true, to a certain extent, with 
respect to discharge current, however, the re- 
sults are not so clearly defined. 

2. The high efficiencies (beam currents) have a 
correspondingly higher percentage of doubly- 
charged ions. 

7 e first of these observations can be explained on the basis of ioni- 
zation cross-sections, while the latter needs further evaluation in terms 
of ion and neutral densities as predicted by beam current and neutral 
loss rate, respectively. 

A linear dependence of the doubly-charged ratio on the ion 
density can be seen in Fig. 20 where neutral loss rate and the 
double ion ratio are both plotted as a function of beam current. For 
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Figure 19. Performance map data relating the double ion ratio to 
discharge mass efficiency. Dashed lines indicate an in 
crease in eV/ion at a given beam current. 
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a thruster operating at maximum propellant utilization the neutral 
loss rate should be a constant regardless of total flow. This, in fact, 
was found to be the case. However, as the neutral loss rate remains 
constant, the double ion fraction increased in a linear fashion up to 
about 2. 0 A beam current* At currents above 2* 0 A the double ion 
ratio began to approach a constant value while the neutral loss rate 
started to increase* It was suspected Lhat at these higher beam cur- 
rents the thruster was operating near the limit of the ion optical sys- 
tem and as a result the mass utilization approached a constant value. 
This, in fact, was shown to be the case by raising the net accelerating 
voltage whereby both the mass utilization and double-ion ratio con- 
tinue' 1 t o increase while the neutral loss rate again approached a con- 
,*lue. Therefore, at a given electron energy and neutral density 
bly-charged ion fraction will increase linearly with ion density. 

If the double -ion production rate is a strong function of 
r|_^ based on the measured thruster values, i. e. , beam current, neu- 
tral loss rate, etc., the same sort of relationships should be observed 

when analyzing the double-ion ratio as a function of thruster radius. 

*f + + 

This case is shown in Fig. 21 where J (r)/J (r) is plotted as a func- 
tion of total ion current density, for beam currents ranging from 1* 0 
to 2. 5 A. It is of interest to note that for any given thruster radial 
location the doubie-ion ratio increases linearly with ion density in 
much the same manner as it does for the thruster with total beam cur- 
rent. However, for constant ion current density the double-ion ratio 
decreases rapidly with increasing distance from the centerline. One 
explanation of this is that Iho neutral density may not be uniform across 
the discharge chamber, but increases with increasing radius. In 
order to further investigate the assumption of a radially increasing 
neutral density, the theoretical analysis of Masek^ is used to predict 
a possible neutral density distribution. The analysis is based on a 
two-dimensional configuration, but since the optics boundary is being 
evaluated, only the solution at z = L. is considered. The analysis 
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Figure 21. Double ion fraction as a function of radial current 
density. 
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predicts an exponential increase in neutral density with radius subject 
to the boundary condition of a known overall neutral loss rate. Once 
the neutral density at the various radial locations is known, it is then 
possible to calculate a radial mass efficiency, n m ( r )> based on the ion 
j\(r), and an equivalent neutral current density J Q ( r )- The results are 
shown in Fig. 22 for the same data points shown in Fig. 21. As can 
be seen the curve is very similar to that of Fig. 19. Thus, although 
the above analysis is approximate, it does indicate that a very nonuni- 
form neutral density distribution could exist, with results consistent 
with overall thruster operating conditions. 

The results of the parametric study can be summarized 
briefly; the double ion ratio was found to be a strong function of the 
mass utilization efficiency, thus, total efficiency, and any gains which 
were made in improving thruster performance were also accompanied 
by an increase in the double ion concentration. This applied in a simi- 
lar manner to the radial parameters and distributions. 

2. Geometric Investigation 

Since the ion and neutral densities are the parameters 
relating mass efficiency and double-charged ion ratio, some assump- 
tions eased on the preceding data can be made with respect to the accel- 
erator grid open area, First, the neutral loss rate varies in direct 
proportion to the accelerator grid open area, and second, if the ion to 
neutral density ratio remains constant, then the double-to- single ion 
ratio will remain constant. The latter assumption is also constrained 
by the fact that the electron temperature and primary-to- maxwellian 
density ratio are not changing appreciably. Based on these assump- 
tions it is possible to define a new mass efficiency 
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Figure 22. Double ion fraction as a function of mass 
efficiency for selected radial locations. 
Data points are. for L = 1.5 A and 
L=2.5A. a 
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where 4? 1 is the new accelerator grid open area fraction and and 

^EM are t ^ ie neutra ^ ^ oss rate an< ^ accelerator grid open area fraction, 
respectively, for the EM optics. In addition to the SHAG optics 
described earlier in this section, a large hole accelerator grid (LHAG) 
set was also assembled to provide additional data on the effects of grid 
open area fractions and thruster performance. The LHAG optics were 
made simply by etching an EM accelerator grid until the desired aper- 
ture size was obtained. This also reduced the thickness, but not 
enough to appreciably affect the structural integrity. The final open 
area fraction achieved was approximately 60%. From eq. (6) an increase 
in mass efficiency should be realized with the use of SHAG optics, pro- 
vided all the assumptions stated earlier are satisfied. Conversely, the 
LHAG were not expected to perform as well relative to EM standards. 

Both sets of optics were evaluated with the results shown 
in Fig. 23. The discharge voltage in all cases was 37 V with discharge 
losses of 185 eV/ion. The mass utilization efficiency was varied by 
changing beam current and mass spectrometer scans were taken at each 
data point. As shown, a 2 to 3% gain in corrected mass efficiency can 
be realized with the use of SHAG optics. However, if one is willing to 
sacrifice part or all of this gain by operation at reduced discharge volt- 
age, then a significant reduction in double ion content can be achieved. 
This is shown in Fig. 24 where the double ion ratio is plotted as a func- 
tion of total, rather than mass, efficiency. By lowering the discharge 
voltage to 3 0 V, a 50% reduction in the doubly-charged ion current can 
be achieved while still maintaining the EM design criteria of 70 to 72% 
total thruster efficiency. Also shown are data points for 33 V discharge 
voltage. 
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Figure 23. Double ion fraction as a function of mass 
efficiency for three different accelerator 
grids. = 37, e T = 185 eV/ion. 










IV* 


THRUSTER POWER CONDITIONER INTERACTIONS 


The interaction between the thruster and the individual power 
supplier of the power processor unit (PPU) needed for its operation 
can significantly affect the capability of either or both to meet the 
exacting specifications that are necessary to make an ion thruster sys- 
tem an attractive alternative to other forms of space propulsion* The 
goal of this portion of the program has been to explore methods for 
minimizing the impact of thruster operating characteristics on the power 
processor requirements that are necessary to establish a stable, reli- 
able system. Additionally, several techniques have been explored for 
reducing the number of power supplies necessary to operate the thruster. 
The following sections describe the work done towards accomplishing 
these goals. 

A. Discharge Current /Voltage Oscillations 

As discussed in Section II, relatively high frequency (10 to 20 kHz) 
oscillations are observed in the discharge voltage or current, dependent 
on the discharge power supply output. The onset and frequency of these 
oscillations is more or less independent of the power supply used, but 
can be related to certain thruster parameters. It was thought initially, 
that it might be possible to find a means of modifying the thruster design 
to eliminate these oscillations. Small variations in several parameters 
produced changes in frequency of the oscillations and/or onset conditions 
but no real progress was made in defining a reliable method of reducing 
or ’eliminating these oscillations. Parameters that were observed to 
influence onset, amplitude, and frequency of discharge voltage /current 
oscillations were: 

• Magnetic baffle coil ampere-turns 

• Cathode propellant flow rate 

• Baffle support geometry 


Baffle diameter 


\ 


All of these parameters suggest that these oscillations are a result of 
some form of current driven instability in the cathode plasma. Two 
types of oscillation were observed. One is characterized by high fre- 
quency oscillation with a relatively stable period and primarily affects 
the discharge current or voltage. The other is a low frequency (50 to 
500 Hz) oscillation that is essentially a relaxation oscillation between 
two discharge modes. The low frequency oscillation affects both the 
beam current and voltage and consequently the entire power processor. 
This low fi'equency oscillation can be damped and essentially elimi- 
nated by increasing the magnetic baffle coil ampere turns. These con- 
ditions were known at the outset of this program and therefore the first 
experiments performed here constituted an attempt to find a physical 
basis for this phenomena. 

Precise measurements of the magnetic field distribution in and 
around the cathode polepiece were made using the precision gaussmeter 
probe positioner shown in Fig. 25* The magnetic baffle current was 
varied and the changes in field distributions were measured and analyzed. 
Figures 26 and 27 show the differences in axial and radial magnetic 
field components for selected locations in thruster S/N301-B (repre- 
sentative of the 700/800 series Engineering Model Thruster), The most 
stable operating points for this thruster are obtained with a magnetic 
baffle current of 3. 5 A ( 10. 5 turn baffle coil). This means that the 
discharge is not subject to "mode shifting" and the high frequency (18 to 
20 kHz) oscillation amplitude is near its minimum value. Note that for 
this baffle current, the axial magnetic field strength goes through zero 
and reverses direction on the axis, inside the pole piece, at a point just 
downstream of the cathode orifice (distance ~3 mm). This measure- 
ment was obtained with a hole in the baffle. Similar axial field mea- 
surements were made using the thruster that had exhibited efficient, 
relatively oscillation free operation during the preceding program 
(NAS 3-16949) and these measurements are shown as Fig. 28* The 
baffle current required to produce this quiescent operation was 4. 5 A 
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Magnetic field plotting apparatus, (a) With alignment fix 
ture. (b) With Gaussmeter probe in place. 
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Figure 26, Variation in axial magnetic field 

component, B z on the thruster axis 
(Thruster S/N301-B) with and without 
magnetic baffle current. 
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Figure 27 . 

Variation in radial magnetic field along 
a radius in the plane of the baffle, with 
and without magnetic baffle current. 
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or greater, and 4, 8 A was the point most frequently run. Note that 
this provides the same kind of axial field condition in the vicinity of 
the cathode orifice as existed in the case for S/N 301 -B, shown in 
Fig. 26* By cutting a slot in the baffle, extensive probe measurements 
were made in the interior of the cathode polepiece in an attempt to 
carefully document the magnetic field conditions thought to be necessary 
for oscillation free operation* Using these field component data, a 
qualitative description was prepared to compare the distribution in 
magnetic field strength and direction in the region between the cathode 
and the main discharge chamber for values of magnetic baffle current 
above and below the "critical” value. This comparison is shown in 
Fig* 29 with the density of lines of force meant to portray field strength* 
Some effort was spent to make this description more quantitative, but 
it was discovered that the slot in the baffle altered the field distribution 
around the annulus region between the baffle and polepiece and perhaps 
in the entire interior* For this thruster, the HRL experimental elec- 
tromagnet thruster, hysteresis was also observed in these magnetic 
field measurements after repeated variations of magnetic baffle cur- 
rent and consequently this approach to investigating thruster discharge 
oscillations was abandoned* Moreover, the cathode propellant flow 
rates for which low oscillation conditions are observed tend to fall 
within a relatively narrow range of values and it is uncertain as to 
which parameter is more important (since they are not independent)* 
Since no really adverse effects could be attributed to the operation of 
a thruster discharge chamber that exhibits discharge cur rent/voltage 
oscillations on the basis of the performance mapping discussed earlier, 
the search for a means of reducing or eliminating this phenomena was 
terminated at this point. 

Here again, as under the previous program, it is concluded that 
several parameters control the transport of electrons within the cathode 
pole region and thereby control the discharge "impedance. 11 Several 
combinations of these parameters may establish electron transport 
to define a discharge "mode" and these parameter combinations may 
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not differ greatly. Thus it is possible to establish conditions so that 
relaxation oscillations between "modes" occur or that an oscillatory- 
behavior is inherent in the electron transport process, or both. To 
eliminate oscillatory behavior consequently is thought to require a truly 
independent means of controlling hollow cathode emission, which may 
or may not be possible, 

B. Beam and Accelerator Power Supply Requirements 

The goal of this project was to determine the feasibility of 
simplifying thruster control by incorporating fast current limiting in 
the beam or screen power supply as an alternative to the high voltage 
turn-off and recycle sequence that was in use at the time of this investi- 
gation. A secondary objective was to evaluate the feasibility of com- 
bining the accelerator power supply with the screen (beam) power supply. 
The investigation began by mapping the voltage-current requirements 
that exist in thruster operation for screen voltages between the off and 
run levels. The experiments performed indicated that current limiting 
alone cannot prevent the thruster from, operating in low voltage, high 
impingement current regions unless the discharge current is reduced. 
The requirement for discharge current reduction to restore high volt- 
age in the screen and accelerator supply is equivalent to a recycle 
sequence and consequently the use of fast current limiting has no real 
advantage over the use of overcurrent initiated high voltage shutoff. 
Moreover, the ion optics assembly sustained damage in the "out of 
perveance" operation that /as permitted in the experiments performed 
to map the voltage-current characteristics and some amount of detri- 
mental operation of this sort can be anticipated with a current limiting 
system. This hazard to the ion optics assembly can be reduced by 
improving the response of the control loops, which also improve the 
dynamic response to load oscillations, etc, A preliminary control 
loop analysis was performed that indicated that a wide bandwidth power 
supply may be practical using lead compensation, but its main advantage 
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would be in dynamic output voltage regulation rather than in 
overcurrent protection by current limiting or high voltage shutoff. 
Investigation of the combined screen and accelerator power supply 
concept was tnore rewarding. 

Many four-trace oscilloscope photographs ofY , V , I , I 

s a s a 

transient events were collected with a variety of HV supply configura- 
tions using the HRL console screen supply and a wide bandwidth pro- 
grammable KEPCO accel supply including: 

1* Screen supply step at 0 accel voltage at various 

discharge currents 

2. Screen supply step with various accel ramp up 
rates at various discharge currents 

3. Accel supply slaved directly to screen supply 
at screen supply step and screen supply ramp 
up at various discharge currents and various 
screen-accel voltage ratios 

4. Accel supply step at 0 screen voltage at various 
discharge currents 

5. Dummy load characteristics of the above 
configurations . 

Conclusions reached from these experiments are: 

1. The screen and accel current, thruster transient 
response is essentially instantaneous. Screen 
and accel current can be defined by the instantane- 
ous screen and accel voltage, 

2, A dual winding HV supply will work with the pres- 
ent thruster design. High voltage application can 
be achieved under a variety of conditions using a 
slaved power supply arrangement, including a V gcrn 
to V acce i ratio of only 100 V to 125 V and from a 
discharge current of 2 A, In all cases the KEPCO 
accel supply current limit value was set for 100 mA, 
realistic for a flight design. In fact, this master- 
slave power supply configuration seems more desir- 
able than independent power supply configurations 
since all intermediate HV operating points between 
HV off and HV fully on were observed to be very 
stable at discharge currents of 2 A or less. The 
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HV could be turned up very slowly with no intermediate 
instabilities or trips. Naturally one would like to 
establish a regime of operating points defining the 
path from thruster off to thruster on where each point 
on the path is unconditionally stable. Unstable regions 
were observed when increasing the discharge current 
after the HV is already established. Engineering 
Model thrusters of the 700/800 series have a tendency 
to oscillate at discharge currents in the range of 
approximately 4 to 6 A with accompanying overload 
trips. Closer inspection revealed that these trips 
result from the relatively low frequency discharge 
oscillations ( —200 Hz) interacting with the screen 
supply output filter causing instantaneous excursions 
of screen voltage which, from conclusion 1, is known 
to be sufficient to drive the thruster out of perveance 
causing a trip. This phenomena is known to be 
improved by additional capacitive filtering in the 
screen supply. The recommendation is that the screen 
supply be respecified for increased filtering. 

3. The backstr earning characteristic of the screen to 

neutralizer current path seems to be (a) monotonic, 

(b) have a definite upper limit of approximately 9 A, 
and (c) behave somewhat like a space charge limited 
diode and finally like a temperature limited diode at 
9 A. 

In summary, the only recommendation that can be made concern- 
ing specification of the screen and accelerator power supply based on 
this investigation is that the accelerator and screen supply can be com- 
bined (or outputs controlled) by a single control, and the frequency 
response of the output voltage regulation should be improved to at 
least 500 Hz . 


C. Utilization of Discharge Power Supply for Auxiliary Functions 

The discharge power supply has the capacity to deliver more 
than 400 W of power and could provide power for at least two auxiliary 
functions that now require separate power supplies. The cathode tip 
and isolator heaters are required only for thruster startup or at very 
low operating power levels. If it is not necessary to operate at low 
power levels for extended time periods, then the discharge power 


supply could be used to supply the preheat function and thereby eliminate 
the requirement for at least two power supplies. The low voltage sec- 
tion of the cathode keeper power supply is another candidate for elim- 
ination* The output of this power supply is not usually varied during 
normal operation, and its function may be satisfied by tapping off the 
desired keeper discharge current with a dropping resistor connected 
to the anoce lead. Again, some flexibility in thruster operation is lost, 
and the power dissipated in the resistor represents a loss, but the 
elimination of the weight and complexity of an extra power supply could 
be sufficient justification to warrant this modification,, 

1, Evaluation of Concepts 

The initial experiments performed under this project 
were planned to demonstrate feasibility of the concept, with minimal 
thruster heater modifications and the emphasis was placed on the cir- 
cuitry for control of the discharge supply. A heater was installed on 
the thruster backplate to simulate the isolator heaters and the dis- 
charge power supply was connected to the thruster as shown in Fig. 30. 
Automatic control algorithms were implemented to operate the heater 
relay, discharge setpoint control and cathode keeper high voltage sup- 
ply as described in Table VII. 

The control algorithm chosen for the relay is very simple — 
being a function only of sensed cathode keeper current. During 
the thruster warmup period the relay is closed and the discharge 
supply is turned on providing power to the series heaters. During this 
period the sensed (in the discharge return leg) and accurately 
describe total heater power. In addition, a clamp circuit is implemented 
in the discharge current setpoint control circuit to automatically limit 
the discharge current to an adjustable predetermined level when the 
relay is closed. This provides a fail-safe means of protecting the 
heater from accidental burn-out and a very simple automatic start up 
algorithm for the discharge supply (i.e. , the discharge supply is always 
on and commanded for the run setpoint). 
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Figure 30. Shared discharge supply wiring diagram. 
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Table VII. Shared Discharge Supply Control Functions 


Function 


Condition 



Cathode keeper HY supply on 


Cathode keeper HV supply off 


CK commanded on and 
< 100 mA 

CK commanded off or 
I^^T* M. > 100 mA 


Discharge control clamped 


IcrT^M, < 100 mA 


Discharge I control unaffected 


> 100 mA 


Ti 873 

At the appropriate time vaporizer power and cathode 
keeper high voltage are turned on (based on automatic control functions 
already existing in the PPU test console). When keeper ignition occurs 
the relay drops out. Presumably, the discharge also ignites and 
arrives at the run setpoint (V^ = 37 V 1^ = 10 A). 

The relay control algorithm was arrived at by consideration 
of all possible states as listed in Table VIII. States 1, 2, 3, 4 
have no significance and state 7 is presumed to be impossible since 
with discharge ignition both low voltage (via the dropping resistor) and 
high voltage are applied to the keeper. For the keeper to fail to ignite 
under these conditions would undoubtedly require a failure in the thruster 
State 5 and state 8 provide the desired states of the relay and require 
no further discussion. State 6, however, requires heater power off in 
which case the cathode keeper will soon starve and return to state 5 
which, presuming proper flow conditions, leads to state 8. In order 
to generate a 1 in state 6 (heater power applied) requires additional 










Table VIII* Heater Relay Logic States 


State 

Discharge on 
C ommand 

Discharge 

Lit 

CK 

lit 

Relay Closed 


X don't care 
X impossible 
X impossible 
X impossible 
1 

if 1 - O. K. htrs. on 
X if 0 - Ck will starve 

and go to state 5 

X impossible 

0 


current sensing in the discharge anode which does not presently exist* 
Otherwise, if one simply bases the decision on sensed I ^ in the return 
teg {where it presently exists) a race will occur in the table between 
state 6 and state 8 because sensed includes heater current* 

Finally, the logic for the cathode keeper HV section is 
shown for illustration only and required no modification, 

A goal of this experiment was to implement the switching 
and control circuitry in a flight qualifiable manner* The system as 
shown impacts the PPU in the following areas: 

1* Addition of cathode keeper dropping resistor 

2. Modification to discharge supply current set- 

point low voltage control logic 

3* Modification of sensed discharge current telem- 

etry to the extent that true discharge current is 
the difference between sensed and sensed 

(note that heater current and discharge plasma 
current are always mutually time exclusive) 
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4. Addition of heater switch element and drive circuitry 

5. Addition of low voltage heater switch control logic. 

It is apparent that 2, 3, and 5 are of no particular concern with respect 
to flight system design* The addition of the dropping resistor is also 
of no particular concern, since for the idea to be valid the resistor must 
have very modest power dissipation (e* g. , 100 ft at 300 mA = 9 W). 
Finally, selection of a switch element to implement the dc power trans- 
fer lead to consideration of two practical choices: (1) a relay, and (2) 

a series pass transistor. Although use of a transistor circuit would, 
in general, have certain difficulties associated with the requirement for 
isolated base drives and the physical mounting considerations associ- 
ated with voltage isolation and power dissipation, for some types of 
missions it may be the best choice. A relay was chosen for expediency 
and as a representative example of a flight qualifiable circuit. In a 
flight system one would expect the stresses on the relay to be rather 
modest in view of the expected number of operating cycles and load 
characteristics of the heaters (mostly resistive). The relay used is a 
Kilovac model HC-1/S43 miniature vacuum relay capable of hot switch- 
ing 18 A and has been used in the 8 cm TYBB program (NAS 3- 17780). 
The power source for the relay coil is assumed to be derivable from 
housekeeping supplies in the PPU and poses no particular design prob- 
lems or penalties. 

Initial experiments were performed by waiting until isolator 
temperature had reached 150°C before initiating vaporizer power. 

The preheat time required for these experiments was on the 
order of 1 hour because the heater current was limited to 5 A by the 
cathode tip heater, and the auxiliary heater was not optimum for heat- 
ing the thruster propellant feed system. Cathode keeper ignition was 
readily obtained using the keeper power supply high voltage section and 
a value of 100 ft for the dropping resistor, R , . Oscilloscope wave- 
forms of the keeper current ignition and steady- state characteristics 
showed excellent, well behaved operation. The keeper was intentionally 
extinguished more than a hundred times followed in each case by suc- 
cessful automatic recycling of heater power and keeper discharge 
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ignition. In addition, many attempts were made to cause loss of keeper 

ignition (while running at = 200 mA from = 150 57 ) by inducing 

various types of high voltage overloads including screen or accel off, 

screen or accel short circuit and screen to accel short from L - 2 A 

13 

steady-state operation. The only loss of keeper and discharge ignition 
observed was for the transient condition with beam (screen) voltage on 
and accel voltage off. Since this is not a normal condition and could not 
be tolerated because of possible ion optics or neutralizer damage, the 
use of the resistor tap on the anode voltage is considered to be a feasi- 
ble alternative to the low voltage power supply now used for steady- state 
keeper operation* Some further experiments were performed to deter- 
mine the limiting value of resistance that can provide this function for 
both start and run conditions. The results showed that: 

1. Reliable, repeatable discharge ignition can be 
forced to occur at vaporizer temperatures 
representing modest flow rates at cathode keeper 
current less than 200 mA, For example, discharge 
ignition occured at rh^ = 80 mA eq, and rh - 1 A 

eq. at 1^,^ “ 105 mA. 

2. Once ignited cathode keeper ignition is maintained 
down to less than 20 mA keeper current over the 
range of test conditions including m^ = — 0 and 

1 m = 0 (T„,, = 242°C and T... , = 243°C) 

m v CV MV ' 

3. The maximum value of the dropping resistor for 
reliable discharge ignition is roughly 150 (I , = 200 mA), 

C iv 

representing about a 6 W power loss, 

4. Reignition of the keeper and discharge after switching 
the discharge and h.v. keeper power supplies off and 
back on again could not be demonstrated with the 1000 
volt high voltage section of the keeper power supply. 
Reignition always required in addition switching on the 
heater and required a relatively random time period 
under the start conditions. 

Thus these preliminary experiments have shown that the 
concept of the switched heater circuit is definitely feasible and will 
likely lead to a much faster warmup sequence and simpler and lighter 
PPU configuration where mission requirements allow mutually exclusive 
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operation of discharge and heaters. To realize the faster warm up 
sequence, a modification of thruster heater will be required so that 
the full discharge power can be utilized. The following section describes 
results obtained with one method of implementing this heater modifica- 
tion. Elimination of the low voltage section of cathode keeper power sup- 
ply is also feasible and thus a high voltage pulsed igniter to supply igni- 
tion also appears more feasible. 

Z. Rapid Start-Up Experiments 

To demonstrate the feasibility of using the discharge 
supply to rapidly heat the thruster it was first necessary to add a new 
heater and modify an existing heater. The new' heater was attached to 
the discharge chamber side of the inner wall of the propellant plenum. 

The size of this heater was selected to match the nominal hot resistance 
of the cathode tip heater. The heater modification consisted of detach- 
ing the end of the cathode isolator heater which is normally connected 
to the thruster shell so that the cathode isolator heater and cathode tip 
heater could be connected in series. This combination was operated in 
parallel with the series combination of the plenum heater and main 
isolator heater as shown in Fig. 31. 

Initially, the test setup was configured as shown in Fig, 30, 
However, while checking the operation of the control logic with the dis- 
charge supply connected to a carbon pile resistor in place of the thrus- 
ter heaters the high voltage relay being used failed. It was then 
derided the basic objective of demonstrating feasibility could be 
achieved by using a separate dc supply switched to replace the discharge 
supply during the heating phase as shov/n in Fig. 31. The power proces- 
sor was also modified to replace the low voltage section of the keeper 
supply with the discharge supply as shown in Fig. 31. 

The operational sequence employed for the rapid start-up tests 
was (1) preheat thruster with commercial dc supply (in lieu of the 
discharge supply), (Z) apply cathode vaporizer and main vaporizer 
power when cathode vaporizer temperature reached approximately 100°C, 

(3) periodically apply cathode keeper voltage until keeper ignition, and 

(4) periodically apply discharge voltage until discharge ignition. The 
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Figure 31. Rapid start-up experiment wiring diagram. 
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power to the isolator, cathode tip, and plenum heaters was removed 
whenever discharge voltage was applied to the anode {see Fig, 31), 

For these tests the output of the dc heater supply was adjusted 
to limit the cathode tip (and cathode isolator) heater current to 4. 5 A. 

At this setting the maximum voltage was 34 V, and 3,8 A was flowing 
through the main Isolator and plenum heaters. The maximum total 
power to the four heaters was 282 W, well below the rating of the dis- 
charge supply for a 30 cm thruster. 

With the sequence described above cathode ignition was achieved 
in 7. 5 to 9 min, and discharge ignition was achieved 11 min after ini- 
tial application of heater power. Temperature profiles at five locations 
on the thruster are presented in Fig. 32 for the first rapid start-up 
test. No effort was made to reduce the start-up time further by optimiz- 
ing the amount of heater power or time of vaporizer power application. 
Keeper operation seemed to be unaffected by the substitution of the dis- 
charge supply for the low voltage section of the keeper supply. 
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THRUSTER- GIMBAL SYSTEM ANALYSIS 

Efforts under this task consisted primarily of (1) a survey of 
gimbal mechanisms vhich have been developed for thrust vector conti ol 
(TVC) under other programs, and (Z) the development of a computer 
program to calculate the,* resultant torques and forces for an arbitrarily 
orientated thruster array* The latter effort was an alternative to the 
original task of designing a specific gimbal actuator and was initiated 
with the concurrence of the NASn Program Manager* 

A. Survey of Gimbal Mechanisms for Thrust Vector Control 

One of the most important factors in thruster array design is 
thrust vector control mechanization. The mechanization approach is 
highly dependent upon requirements such as thrust direction accuracy, 
coupling between the three torque axes, thrust loss, thruster particle 
(ions and atoms) flux distributions, reliability and the particular mis- 
sion under consideration. For example, for planetary missions, cruise 
mode perturbations and biases are relatively small and thruster gimbal 
angles of 5 to 10° provide sufficient control authority* However, 
on out-bound type trajectories, the number of operating thrusters is 
gradually decreased. In order to maintain the thrust vector through 
the center-of -mass, additional gimbaling is required. Thus, the gim- 
bal requirement for a planetary mission design must include the vari- 
able number of operating thrusters as well as control authority. 

In the case of earth orbital (EO) missions, the number of operat- 
ing thrusters is generally less variable since available power is rela- 
tively constant. (An exception to this condition occurs if orbit raising 
through the Van Allen Belt is considered. ) However, the torque required 
for EO mission maneuvers {e, g. , vehicle attitude changes of 180° twice 
per orbit) is one to two orders of magnitude greater than for planetary 
missions. The gimbal angles needed to produce such torcues are gen- 
erally in the range of 20 to 30° and thus result in significant thrust 
losse s . 
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The first TYC hardware developed was designed to minimize 
as many uncertainties as possible by use of a two axis array translator 
and single axis thruster gimbals. This approach provided relatively 
accurate knowledge of thrust vector pointing, good control authority, 
and minimized coupling between axes. Unfortunately, the translator 
single point failure possibilities, combined with its weight and com- 
plexity resulted in a search for alternatives* 

Once the thrust vector accuracy requirements were relaxed 
and the elimination of single point failures were encouraged, the logical 
choice was individual thruster two axis gimbaling. 

Workers at Jet Propulsion Laboratories (JPL) carried this 
design concept into a prototype hardware phase* The JPL concept, 
illustrated in Fig, 33, packages both gimbal actuators into a single 
housing (Fig. 34). To provide for launch vibration, a pyro- released 
latch is used to support the thruster free gimbal nad. The JPL overall 
packaging approach was selected to provide actuator thermal control, 
to prevent thruster heat rejection to the vehicle, and for thruster ther- 
mal control. Based on analytical studies, the open sides and slanted 
thermal blanket were believed to be necessary for thruster thermal 
control. The JPL "thrust assembly” designed for the Encke Flyby 
Mission^ was estimated to weigh 91 kg, or 15.2 kg /thruster. 

Recent studies of the Solar Electric Propulsion State (SEPS) by 
13 14 

Rockwell International and Boeing for Marshall Space Flight Center, 
included thruster array conceptual design for up to twelve thrusters. 

For these studies, thruster array packaging constraints, thermal con- 
trol limits, and interface definition were provided by Hughes Research 
Laboratories and NASA LeRC Research Center, The resulting designs 
were relatively strongly influenced by the results of tests performed by 
Mirtich and Rawlin of LeRC. In the LeRC tests, it was concluded that 
thrusters could be stored at temperatures down to -200°C, and that 
thruster operation was essentially unaffected by thermally insulated 
surrounding or by direct solar radiation up to 2.5 sun 1 s intensity , Thus, 
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Figure 33. JPL, gimbaling concept. 
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these designs reflect approaches in which thruster thermal control is 
not required. However, both approaches provide insulation between 
the thrusters and vehicle body. 

Both Rockwell and Boeing selected thruster two axis gimbaling 
concepts. The Rockwell concept, similar to that shown in Fig. 34 uses 
linear actuators of the type developed at JPL for midcourse rocket 
engine gimbaling. The Boeing concept, shown in Fig. 35, uses motor 
driven tape type actuators to rotate the thruster around the thrust axis 
and around the thruster gimbal axis. The tape drive approach poten- 
tially provides simplicity and low backlash, but may require shielding 
from direct solar radiation, 

A comparison of the failure rates and weights for these three 
systems as well as the more conventional tape driven mechanism 
(Fig. 36) has been made by breaking the systems up into similar com- 
ponents as shown in Table IX. By using consistent failure rates for 
similar components, the comparison should be valid even though each 
system uses different gimbaling techniques. On the basis of the simi- 
larity in failure rates seen here, selection of a gimbal mechanism on 
the basis of reliability is not considered possible, however the Lower 
weights for the Boeing design and the conventional tape drive mecha- 
nism favors the latter two designs. The weights shown in Table IX 
are for a single thruster and were estimated using consistent weights 
for similar components. 

B* Thruster Gimbaling Computer Program 

A necessary step in the general area of gimbal requirement 
definition and gimbal approach selection is the generation of para- 
metric trade study data. In order to select a given thruster array pat- 
tern and gimbal angle, one must have the means for computing torques 
about three axes, the maximum or net gimbal, and thrust losses. In 
addition, the gimbal angle needed to achieve a given torque level, will 
place constaints on thruster-to-thruste r positioning. Thruster position- 
ing in *;urn will affect torque capability, making the selection of a gim- 
bal angle and thruster pattern an iterative process. To facilitate this 
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Figure 35. Boeing tape drive gimbaling concept. 
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Figure 36. Conventional tape drive gimbaling concept 




Table IX. Summary of Estimated Failure Rates (failures/ 10 hr) 
and Weights for Different Thruster Gimbal Systems 


JPL 

Boeing 

Rockwell 

Component Two Axis 

Tape 

Linear 

Activation 

Drive 

Actuator 


Stepper Motors 0.22 0. 2Z 

Bearings, Main 0. 1 1 0. 1 1 

Shaft, Stub 0. 16 0. 08 

Shaft, Main Moi nt — 0, 14 

Potentiometers 0.25 0.25 

Metal Tapes 0. 18 0. 18 

Pulleys 0. 14 0. 14 

Ball Screws — — 

Rod Ends, Bearing — — 

Bearings, Motor Gimbal — — 

Gear Sets 0,018 0.018 

Bearings, Gear 0.008 0.008 

Shaft Seal — — 

Gasket — — 

Bell Cranks — — 

Leaf Spring Links — — 

Total Failure Rate 1.09 1. 15 

Estimated weight (Kg) 3.26 2.35 


0 . 22 
0 . 11 

0. 28 
0. 25 


0. 07 


0. 018 
0. 008 
0. 03 
0 . 01 

0. 17 

0 . 20 

1. 37 
3.9 



0 . 22 
0 . 1 1 
0. 48 

0. 25 


0. 07 
0. 02 
0 . 22 
0. 018 
0. 008 


1. 40 
2. 4 


process a computer program was written to calculate thrust, total 
gimbal angle for each thruster, net thrust including gimbaling, thrust 
vector direction, and torques about three orthogonal axes as a function 
of the following parameters: 

Number of operating thrusters 

Thruster pattern 

Vehicle center-of-mass location 









Gimbal angles (or and p ) for each thruster 

Gimbaling strategy (i.e., choice of which thrusters 
are to be used for roll, pitch, and yaw functions). 

The computer program outputs are: 

Total thrust loss 

Resultant thrust vector 

Resultant torques 

Net gimbaling angle for each thruster, 

C. Analysis 

!♦ Forces 

In vector rotation, the thrust vector of the ith thruster 
in the vehicle coordinate system is related to the thrust vector in the 
primed (thruster) coordinate system by - \F/ where 



and \ is the matrix of direction cosines which represent two successive 

rotations by a and (3, For the case shown in Fig. 37, assuming rotation 

by angle (3 about this Y axis is applied first, X. is given by the product 

of matrices representing rotations about the x and y axis. Thus ior a 

15 

rotation by angle (3 about the Y axis the matrix is 


B = 


1 G 0 
0 co s p sinp 
0 -sin(3 cos(3 


( 8 ) 







and for a rotation by angle $ about the x axis the matrix 


A = 


( coso 1 sin ca 
-sine* coso 1 
0 0 



( 9 ) 


The resultant matrix for two successive rotations assuming 
(3 is applied first is given by the product of these two matrices or 


= AB - 


( 


cosa (sintf cosj3) 
-sinor (cose* cosp) 
0 - sinp 


sina 

COScJ) 

cos 



( 10 ) 


Since there is only a Z component of force in the primed 
system eq* (7) can be written in component form in the vehicle coordi- 
nate system using eqs. (7) and { 10) (for each thruster) 

£ - (cosa) (cos (3) I f I 

f - (sina) (cos (3) | f | (11) 

£ z = ( - sinp) | f | 


The x and y gimbal angles a and (3 combine to form a total gimbal angle 
Y given by 


l f 


= cos 


(cosa cosp) 


( 12 ) 


where [ f j is the magnitude of the total force for each thruster. The 

resultant total force for the thruster array is evaluated by summing 

eq. (11) for each thruster. That is, the total forces F , F , and F 
^ x y z 


given by 


are 


) 

i 


f = £f . 

X * 7 * X 
1 


F y = ? f y • 


F = £f 

Z ~ z 
1 


The resultant thrust loss is 


Thrust loss - 


1 " sT / 100 * (%) 


The resultant thrust vector makes angles A, B, and C 
with the x,y, and z axes of the vehicle coordinate system and can be 
calculated as follows: 


A = cos ^ ~ 

f t 


B = cos 


l_v 


where 


C = cos 


' ? 7 7 

, / F + F + F 
y x y z 
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Torque s 


In vector notation the torque about the vehicle CG due 
to an arbitrary thruster is given by (see Fig* 35) 

T = RO x F (17) 


and 


RO - RT - RG . (18) 

So in component form, if a thruster has position components x 1 , yl, 
and zl, in the vehicle coordinate system and the CG has coordinates 
xo, yo, and zo, then using (3 8) the component for the RO vector are 
given by 


r = xl - xo 
x 


r = y 1 - yo 

y 7 


(19) 


r = z 1 - zo 
z 


So in component form the torques from an individual thruster are: 


t - r f - r f 
x y z z y 


t = r f - r f 
y z x x z 


( 20 ) 


t “ r f - r f 
z x y y x 


where f , f , and f z are given by eq, (11)* The resultant total torques 
x y 

are evaluated by summing eq. (2 0) for the array. 
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3. Computer Program 


A Fortran IV computer program was written to solve 
the above torque and force equations. The program was developed on 
the DEC-system 10 time sharing system and should run on other com- 
puter systems with minor coding changes. A complete program source 
listing is shown in Appendix A, To illustrate the use of this program we 
will solve for the resultant torques, forces etc., for - 3 x 3 thruster 
array with seven operating thrusters. For this configuration the two 
corner pairs and an inner pair are gimbaled to produce the maximum 
roll torque subject to the constraint of a maximum total gimbal angle 
of 28° for any gimbalied thruster. These thrusters are mounted on a 
27 in. (0* 636 m) square pattern and the vehicle CG is symmetrically 
placed 56 in. (1,422 m) to the rear of the thrusters mounting plate. The 
computer print-out for this configuration is shown in Fig. 38. The program 
first prints out the input data which include the thruster and CG coordi- 
nate and the gimbal angles. The output data includes the thrust and 
torque components for each thruster along with the total gimbal angles 
and a summary of resultant torques and forces and the total thrust loss, 

7 THHUSTliriS* Co AT X= 0.00M Y = O.GOT Z=-!*42 


THHUSTtitf XUO Y(M) ZU) ALP(OhG) BET(DEG) TrlRUST(KLBS) 


1 

0 . 69 

0.69 

0. 00 

20.0 

20.0 


30.0 



-0.69 

0.69 

0.00 

-20.0 

23 . 0 


30.0 


3 

-0.69 

0.00 

0.00 

23.0 

0.0 


30.0 

h \ 

4 

0.00 

0.00 

0.00 

0.0 

0.0 


30.0 


5 

0.69 

0.00 

0 . 00 

28.0 

0.0 


30.0 


6 

0.69 

-0.69 

0.00 

20.0 

-20.0 


30.0 


7 

-0.69 

-0.69 

0.00 

-20.0 

-20.0 


30.0 


THRUSTER 

FX 

FY 

FZ 

TX 

'l'Y 

TZ 

THETA 


1 

-10.26 

9.64 

26.49 

.015 

-.107 

.045 

27.99 


2 

-10.26 

-9.64 

26.49 

. 106 

.012 

.045 

27. 99 


3 

0.00 

14.08 

26.49 

-.066 

.060 

-.032 

28.00 


4 

0.00 

0.00 

30.00 

.000 

. 000 

.000 

0.00 


5 

0.00 

14.08 

26.49 

-.066 

-.060 

.032 

28 . 00 


6 

10.26 

9.64 

26.49 

-. 105 

-.012 

.045 

27.99 


7 

10.26 

-9.64 

26.49 

-.015 

. 107 

.045 

2 7.99 


TOTAL: 

0. 00 

28.17 

l 33 - 94 

-0.131 

-0.000 

0. 1 / 9 



3TUUST LOSS- 

10.03% 

A= 90.00deq L3= 

8 l . 52deq 

U= 8. 

48deq 


: 3 


Figure 38. Computer printout for thruster gimbaling computer program. 
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VI. CONCLUSIONS 

The work performed under this 2. 5 kW Advanced Technology Ion 
Thruster Program provided a quantitative assessment of the perfor- 
mance capabilities and limitations of the 700/800 series Engineering 
Model Ion Thruster, Measurements made over a range of thruster 
control parameters show that the total thruster efficiency is not changed 
significantly by parameter selection if the data is corrected to account 
for thrust losses from doubly-charged ions and beam divergence. Per- 
formance data were obtained for variations of beam current (throttling) 
from 0. 5 A to 2. 75 A and the corrected thruster efficiency is in good 
agreement with that predicted using Kaufman's maximum propellant 
utilization theory. Equilibrium temperatures of thruster components 
were well within the limits considered satisfactory except at the lowest 
beam current levels where r,ome isolator heater power is required to 
prevent propellant condensation. 

Double ionization is a major problem area and places limitations 
on both the thruster efficiency and the discharge chamber lifetime. 

Data obtained under this program show a definite relationship between 
doubly-charged ion production and propellant utilization efficiency (for 
a given thruster). The transmission of the ion optics accelerator grid 
is found to have the greatest effect on this proportionality and reduction 
in accelerator aperture diameter was the only method observed to be 
successful in significantly reducing double ionization. The capabilities 
ana/or the disadvantages of this design modification were by no means 
exhaustively explored under this program, and many questions remain 
to be answered* Discharge voltage, discharge power and beam current 
also influence doubly-charged ion production, but to lessor degrees. 
High power operation obtai led by increasing beam current using this 
thruster design will definitely reduce thruster lifetime because of 
increased double ionization. 
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Thruster /power conditioning interactions were explored and 
concluded to be more a question of the specific thruster and power 
supply combination than a general problem. Oscillations in discharge 
current and/or voltage are nearly always observed. These oscillations 
were not proven to be detrimental in any manner that affects thruster 
vf r L;iency or overall thruster operability. Isolated incidents of trouble- 
some operation have been observed, but in general, the power condi- 
tioning as defined under preceding programs (Contracts NAS 3-14140, 
NAS 3- 16528, NAS 3- 16949) can provide satisfactory operation and con- 
trol for the 700-800 series Engineering Model Thruster ove-r a wide 
operating range. Several options for using the discharge power supply 
to provide the functions of other separate power supplies were demon- 
strated and consequently the power conditioning could be made lighter 
and more reliable by eliminating these separate power supplies. Simi- 
larly, t was demonstrated that the accelerator supply could be incor- 
porated m the screen (beam) power supply, 

A study of gimbal actuator requirements showed that a ‘'universal 
actuator design is not feasible because of the relatively wide range of 
requirements that can be identified with specific applications. A com- 
puter program was developed to aid in analyzing gimbal actuator require 
merits in place of attempting to design a general purpose actuator. 


102 


REFERENCES 


1. R. L. Poeschel, H. J. King, and D. E. Schnelker, "An 
Engineering Model 30 CM Thruster, " AIAA Paper No. 73-1084, 
Lake Tahoe, Nev. (1973). 

2. J. S. Sovey and H. J. King, "Status of 30 CM Mercury Ion 
Thruster Development" AIAA Paper No. 74-1117, San Diego, 

Calif. (1974). 

3. H. J. King, et al. , "Low Voltage 30 CM Ion Thruster Develop- 
ment, " NASA CR- 13473 1, October 1974. 

4. B. G. Herron, J. D. Worden, and J. M. Simpkins, "A 30 CM 
Thruster Power Processor Test Console, 11 AIAA Paper No. 

73- 1104, October 1973. 

5. B. Banks, V. Rawlin, A. Weigand, and J. Walker , n Direct 
Thrust Measurement of a 30 cm Ion Thruster." 

6. H. R. Kaufman, "Ion- Thruster Propellant Utilization, " 

J. Spacecraft and Rockets _9, No. 7, July 1972, pp. 511-517. 

7. F. F. Terdan and R. T. Bechtel, "Control of a 30 CM Diameter 
Mercury Bombardment Thruster, " AIAA Paper No, 73-1079, 

Lake Tahoe, Nev. (1973). 

8. HRL Technical Staff, "2. 5 kW Advanced Technology Ion 
Thruster, " NASA CR- 134687, August 1974. 

9. R. P. Vahrenkamp, "Measurement of Double Charged Ions in 
the Beam of a 30 cm Mercury Bombardment Thruster, " AIAA 
Paper 73- 1057, Lake Tahoe, Nevada, October 1973. 

10. C. Collett, "A 7700 Hour Endurance Test of a 30 CM Kaufman 
Thruster, " AIAA Paper No. 75-366, New Orleans, La. , 

March 1975. 

11. T. D. Masek, "Plasma Properties and Performance of 
Mercury Ion Thrusters," AIAA Paper 69-256, March 1969. 

12. G. S. Perkins, et al. , "A Mechanism for Three-Axis Control 
of an Ion Thruster Array, " AIAA Paper 70-1156, AIAA 8th 
Electric Propulsion Conference, Stanford, California, August, 
1970. 

13. J. H. Duxbury, and J. H. Gerpheide, "Solar Electric Propul- 
sion Encke Slow -Flyby 1979 Mission and Spacecraft Description, " 
JPL Report 701-200, June 28, 1974. 


103 


"Concept Definition and. Systems Analysis Study for a Solar 
Electric Propulsion Stage, " Rockwell International Final 
Report on NAS 8-30920, SD 74-SA-0176, February 1975. 

"Concept Definition and Systems Analysis Study for a Solar 
Electric Propulsion Stage, " Boeing Aerospace Company 
Final Report on NAS 8-30921, D180-18553, January, 1975. 

H. Goldstein, "Classical Mechanics," Addison and Wesley, 
1959, p. 126. 


APPENDIX 


SOURCE LISTING FOR THRUSTER GIMBALLING 
COMPUTER PROGRAM 


A source program listing for the thruster gimballing computer 
program is shown below. The input data for the thruster and CG coor 
dinates are in meters, the forces are in mlbs and the gimbal angles 
are in degrees. The program outputs are given in Ib-ft for the torque 
and mlbs for the forces. 


DIMENSION XX(ZO), YY(20), ZZ(20), AA(20), BB(20), CC(20), 
TH(20 )OPEN(UNIT=2 0, FILE= 1 DATA. DAT') 

C READ # OF CASES, # OF THRUSTERS, AND CG COORDINATES 
(IN METERS) 

READ(20,*) NCASE, NTHRUS, X), YO, ZO 
WRITE(5, 100) NTHRUS, XO, YO, ZO 
C 1 = 3. 14159/180. 

C2 = . 00328084 
DO 20 JXX= 1, NCASE 
C C2 CONVERTS MLBS*M TO LB*FT 

WRITE(5, 2 00) 

DO 10 1=1, NTHRUS 

C READ THRUSTER COORD. , GIMBAL ANGLES, AND THRUST 
(M, DEG, AND MLBS) 

READ(20, *) XX(I), YY(I), ZZ(I), AA(I), BB(I), TH(I) 

10 WRITE( 5 , 3 00) ( , XX(I), YY(I), ZZ(I), AA(I), BB(I), TH(I) 

C INITIALIZE SUMS 

TT = 0. 

FXX-0. 

FYY=0. 

FZZ = 0. 

TXX = 0. 

TYY=0. 

TZZ = 0. 

WRITE(5, 310) 

WRITE(5, 550) 

DO 6 1=1, NTHRUS 
X=XX(I) 

Y=YY(I) 

Z = ZZ(I) 

Al= AA(I) 

B 1=BB(I) 

T = TH(I) 

tt_tt+t 


CAL GIMBAL ANGLES IN RADIANS 
A=C1*A1 
B=C1*B1 

COMPUTE COMPONENTS OF VECTOR BETWEEN CG AND 
THRUSTER 

RX=(X-XO)*C2 
RY=( Y=YO)*C2 
RZ = (Z-ZO)*C2 

COMPUTE FORCE COMPONENTS IN S/C COORD. SYSTEM 
FX = -SIN(B)*T 
FY=SIN(A)*COS(B)*T 
FZ=COS(A)*COS(B }*T 
COMPUTE TORQUES ABOUT CG 
TX=RY*FZ-RZ*FY 
TY=RZ*FX-RX*FZ 
TZ=RX*F Y-RY^FX 
ADD UP TORQUES AND FORCES 
TXX=TXX+TX 
TYY=TYY+T Y 
TZZ=TZZ+TZ 
FXX=FXX+FX 
FYY=FYY+FY 
FZ Z = FZZ+FZ 

CALCULATE TOTAL GIMBAL ANGLE FOR EACH THRUSTER 
(IN DEG) 

TOT=ACOS(COS(A)*COS(B))/C 1 
6 WRITE(5, 400) I, FX, FY, FZ , TX, TY, TZ, TOT 
FTOT=SQRT)FXX + FXX+FYY*FYY+FZZ*FZZ) 
THETAX=ACOS(FXX /FTOT)/CI 
THETAZ = ACOS(FZZ /FTOT)/C 1 
THETA Y = ACOS(FYY /FTOT )/C 1 
TLOSS=( 1 . -FZZ/TT)*100. 

WRITE(5, 310) 

WRITE{5, 500) FXX, FYY, FZZ, TXX, TYY, TZZ 
WRITE(5, 350) TLOSS, THETAX, THETAY, THETAZ 
20 CONTINUE 

500 FORMA T(6X, 'TOTAL: ', F6.2,2X, F6.2, IS, F6.2 & 

2X, F6. 3, 3X, F6. 3, 2X, F6. 3) 

500 FORMAT(6X, 'THRUSTER FX FY FZ 
& .TX TY TZ THETA 1 ) 

100 FORMAT( 12X, 13, 'THRUSTERS: CG AT X=\F5.2, 'M Y=', 

& F5.2, 'M z = ',F5.2) 

200 FORMAT(lH-, 6X, 'THRUSTER (X(M) Y(M) Z(M) ALP(DEG) 
BET(DEG) &: THRUST (MLBS)') 

310 FORMAT) 1 5X, ' ') 

300 FORMA T(6X, 13, 6X, F5. 2, 3X, F5.2,2X, F5.2, 3X, F5. 1, 
4X,F5. 1, 8X,F5. 1) 

400 FORMAT) 6X, 13, 5X, F6.2,2X, F6. 2,2X, F5. 2, 3X, F5. 3, 

4X, F5. 3, 3X, &F5. 3, 3X,F6.2) 

350 FORMAT) IX, 'THRUST LOSS=',F6.2, '% A=',F6.2, 'deg 
& B = \ F6.2, 'deg G=',F6.2, 'deg') 

END 


